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Abstract: Paragonimiasis is an important food-borne parasitic zoonosis caused by infection with lung flukes of the genus 
Paragonimus. Of the 7 members of the genus known in Thailand until recently, only P. heterotremus has been confirmed 
as causing human disease. An 8th species, P. pseudoheterotremus, has recently been proposed from Thailand, and has 
been found in humans. Molecular data place this species as a sister species to P. heterotremus, and it is likely that P. pseu- 
doheterotremus is not specifically distinct from P. heterotremus. In this study, we collected metacercariae of both nominal 
species (identification based on metacercarial morphology) from freshwater crabs from Phetchabun Province in northern 
Thailand, Saraburi Province in central Thailand, and Surat Thani Province in southern Thailand. In addition, we purchased 
freshwater crabs imported from Myanmar at Myawaddy Province, western Thailand, close to the Myanmar-Thailand bor- 
der. The DNAs extracted from excysted metacercariae were PCR-amplified and sequenced for ITS2 and cox7 genes. The 
ITS2 sequences were nearly identical among all samples (99-1 00%). Phylogenies inferred from all available partial coxl 
sequences contained several clusters. Sequences from Indian P. heterotremus formed a sister group to sequences from 
P. pseudoheterotremus-type metacercariae. Sequences of P. heterotremus from Thailand, Vietnam, and China formed a 
separate distinct clade. One metacercaria from Phitsanulok Province was distinct from all others. There is clearly consid- 
erable genetic variation in the P. heterotremus complex in Thailand and the form referred to as P. pseudoheterotremus is 
widely distributed in Thailand and the Thai-Myanmar border region. 
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INTRODUCTION 

Paragonimiasis is an important food-borne zoonosis, espe- 
cially in Asian countries, caused by infection with lung flukes 
of the genus Paragonimus. More than 50 species of Paragonimus 
have been described from Asia, America, and Africa [1]. Seven 
species, namely, P. heterotremus, P. westermani, P. siamensis, P. 
bangkokensis, P. harinasutai, P. paishuihoensis, and P. macrorchis - 
have been reported from Thailand. Of these, only P. heterotre- 
mus was regarded until recentiy as causing human disease in 
Thailand [2,3], as it also does in India [4-6], Vietnam [7-9], 
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and Laos [10]. Recently, an 8th nominal species, P. pseudohet- 
erotremus, was described from Thailand [11], and human pul- 
monary paragonimiasis caused by this species has been identi- 
fied using DNA sequence data of the internal transcribed spac- 
er 2 (ITS2) and cytochrome c oxidase subunit 1 (coxl) from 
eggs in the sputum of a patient [12]. This species was distin- 
guished from P. heterotremus by the small size of its metacer- 
cariae and by details of adult surface spination. P. pseudohetero- 
tremus appears, on the basis of DNA sequence data, to be a 
very close sister species to P. heterotremus. In Vietnam, metacer- 
cariae resembling those of both of these species have been 
found, but cannot be separated using molecular data (Doanh, 
pers. comm.). P. pseudoheterotremus may thus not be specifical- 
ly distinct from P. heterotremus. However, both names will be 
used in this paper, simply to indicate the size differences in 
metacercariae. 

The discovery of a human case attributed to P. pseudohetero- 
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tremus [12] prompted us to study the distributions and rela- 
tionships of metacercariae of P. heterotremus and P. pseudohet- 
erotremus types in freshwater crabs in Thailand. Moreover, nei- 
ther of these types were reported from Myanmar, which shares 
a long border with Thailand. Here, we report ITS2 and coxl gene 
sequences of metacercariae of P. heterotremus and P. pseudohet- 
erotremus types from several regions in Thailand and an area in 
Myanmar, and place these in a phylogenetic context using new 
and published sequences. 

MATERIALS AND METHODS 

Collection of Paragonimus metacercariae 

Paragonimus metacercariae (Fig. 1) were collected from natu- 
rally infected freshwater crabs which were collected in 4 prov- 
inces of Thailand (Fig. 2); Phitsanulok Province and Phetch- 
abun Province in the north, Saraburi Province in the center, 
and Surat Thani Province in the south of Thailand. In addi- 
tion, crabs were also collected in Myawaddy Province, Myan- 
mar, bordering western Thailand. The crabs were digested with 
artificial gastric juice (1% pepsin in HC1 solution). Digested 
samples were filtered through wire sieves (pore size; 300 and 
500 pm), and the filtrate was allowed to stand for sedimenta- 
tion. The sediments were examined under a stereomicroscope. 
The metacercariae found were excysted using a method previ- 
ously described [13]. Each freshly excysted Paragonimus meta- 
cercaria was morphologically identified using microscopy be- 
fore DNA extraction. 





Fig. 2. Locations of provinces in Thailand and Myanmar from 
which freshwater crabs were sampled for this study. Myawaddy: 
Myanmar (1 ; present study), Thailand: Phitsanulok (2; present 
study), Phetchabun (3; present study), Saraburi (4; present study 
and Thaenkham and Waikagul, 2008), and Surat Thani (5; present 
study). In addition, provinces from which P. pseudoheterotremus 
has been reported in previous studies are shown: Kanchanaburi (6; 
Thaenkham and Waikagul, 2008) and Loei (7; Intapan et al., 2012). 




Fig. 1 . Representative metacercariae morphologically resembling those of P. heterotremus (A) and P. pseudoheterotremus (B). The P. 
heterotremus-type metacercaria (A) was collected in Phisanulok Province, northern Thailand, whereas the P. pseudoheterotremus-type 
metacercariae (B) were collected in Surat Thani Province, southern Thailand. 
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DNA extraction 

To extract genomic DNA, each excysted Paragonimus meta- 
cercaria was homogenized separately with a disposable poly- 
propylene pestle (Bellco Glass Inc., Vineland, New Jersy, USA). 
Then, DNA was extracted using a NucleoSpin Tissue Kit (Ma- 
cherey-Nagel GmbH & Co., Diiren, Germany). The DNA was 
eluted in 100 pi of distilled water and stored at -20°C until used. 

PCR, sequencing, and molecular phylogenetic analysis 

The ITS2 sequence was amplified using primers PhITS2-F 
(5'-CTG TGT GAA TTA ATG TGA ACT GC-3') and PhITS2-R 
(5'-AGT GAT ATG CTT AAG TTC AGC G-3') [12], which were 
designed from the known 1TS2 sequence of P. heterotremus from 
northeastern India (GenBank accession no. AB308377). A 
fragment of the coxl gene was amplified using primers, Pph- 
COI-F (5-CCG GGTTTG GTG TTG TG-3') and PphCOI-R (5'- 



Table 1 . Accession numbers for new Paragonimus sequences 
deposited in GenBank 



Country 
(Province) 


cox 1 


ITS2 


Sequence 
ID 


GenBank 
No. 


Sequence 
ID 


GenBank 
No. 


Myanmar: 


PhC01MW4 


KC859926 


Phl2MW4 


KC894640 


(Myawaddy, 


PphC01MW5 KC859936 


Pphl2MW5 


KC894649 


MW) 


PphC01MW6 KC859937 


Pphl2MW6 


KC894650 


Thailand: (Phit- 


PhC01PL2 


KC859931 


PN2PL2 


KC894644 


sanulok, PL) 










Thailand: 


PhC01PB3 


KC859927 


PN2PB3 


KC894641 


(Phetchabun, 
PB) 










Thailand: 


PhC01SB1 


KC859933 


Phl2SB1 


KC894646 


(Saraburi, SB) 












PhC01SB3 


KC859935 


Phl2SB3 


KC894648 


Thailand: (Surat 


PphC01ST1 


KC859938 


Pphl2ST1 


KC894651 


Thani, ST) 


PphC01ST2 


KC859939 


Pphl2ST2 


KC894652 




PphC01ST3 


KC859940 


Pphl2ST3 


KC894653 




PphC01ST4 


KC859941 


Pphl2ST4 


KC894654 




PphC01ST5 


KC859942 


Pphl2ST5 


KC894655 




PphC01ST6 


KC859943 


Pphl2ST6 


KC894656 




PphC01ST7 


KC859944 


Pphl2ST7 


KC894657 




PphC01ST8 


KC859945 


Pphl2ST8 


KC894658 




PphC01ST9 


KC859946 


Pphl2ST9 


KC894659 



ACA ACG AAC CAA GTG TCA TG-3') [12], which were de- 
signed from a region of P. pseudoheterotremus coxl gene (Gen- 
Bank accession no. EF446315). 

PCR was conducted using a GeneAmp® PCR System 9700 
(Applied Biosystems, Singapore). The reaction was carried out 
in a 25 pi volume containing PCR 1 x FastStart High Fidelity 
Reaction Buffer (Roche Applied Science, Mannheim, Germa- 
ny), 1.8 mM MgCh, 0.2 mM of each deoxyribonucleotide tti- 
phosphate, 0.2 pM of each primer, and 0.625 units of FastStart 
High Fidelity Enzyme Blend (Roche Applied Science). The DNA 
template was initially denatured at 94°C for 5 min. The ampli- 
fication procedure comprised 35 cycles at 95°C for 30 sec (de- 
naturation), 55°C for 30 sec (annealing), and 72°C for 30 sec 
(extension), with a final extension at 72°C for 10 min. The am- 
plified product was subjected to electrophoresis on a 1.5% 
agarose gel. The PCR product was then cut and purified for 
DNA sequencing, which was performed using the MegaBACE 
1,000 DNA Analysis System (GE Healthcare, Piscataway, New 
Jersy, USA). 

The partial ITS2 and coxl gene sequences of individual Para- 
gonimus metacercariae from each locality were deposited in 
GenBank under accession numbers presented in Table 1. They 
were analyzed using the BLAST-N search (National Center for 
Biotechnology Information, Bethesda, Maryland, USA). Pub- 
lished coxl sequences from P. heterotremus and P. pseudohetero- 
tremus were aligned with our new sequences (alignment length 
309 bp when uimmed to the length of the shortest sequence) 
using ClustalW [14] and a maximum likelihood tree construct- 
ed using MEGA 5.2 [15]. The best-fit substitution model was 
determined in MEGA to be the Hasegawa-Kishino-Yano (HKY) 
model with uniform rates among sites, but assuming a pro- 
portion (0.65) of invariant sites. Sequences from P. westermani 
were used to provide an outgroup. 

RESULTS 

The partial ITS2 sequences of all metacercariae from differ- 



Table 2. Pairwise percentage differences in range (mean) between clusters on the tree 





(2) 


(3) 


(4) 


(5) 


P. heterotremus (excl. India) (1) 


7.8-9.1 (8.1) 


8.7-11 (9.5) 


7.8-9.1 (8.3) 


13-20 (16.3) 


P. heterotremus (India) (2) 




4.5-6.5 (5.5) 


8.1-9.1 (8.5) 


14-20 (16.8) 


P. pseudoheterotremus (3) 






9.1-10(9.9) 


14-21 (16.7) 


P. heterotremus Phitsanulok (4) 








16-21 (18.1) 


P. westermani (5) 
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100 
70 



18 



65 



PhC01SB1 KC859933 

P. heterotremus AB354230 THA 
- P. heterotremus EF446319 THA 
P. heterotremus EF446321THA 



P. heterotremus EF446318 THA 
PhC01SB3 KC859935 
P. heterotremus HM627189 CHN 
P. heterotremus AB354229 THA 
P heterotremus AY618839 CHN 
. P. heterotremus AF159597 THA 
"I P. heterotremus AY618840 CHN 
P. heterotremus AB270675 VNM 
P. heterotremus AY6 18842 VNM 
P. heterotremus AY618843 VNM 
I P. heterotremus EF446320 THA 
_ PhCCM MW4 KC859926 
, P. heterotremus AB270676 VNM 
I P heterotremus AB270677 VNM 
PhC01PB3 KC859927 
. P. heterotremus HM627190 CHN 
P. heterotremus HM627188 CHN 
P. heterotremus AY618841 VNM 
P. heterotremus DQ234301 VNM 
-PhC01PL2 KC859931 

P. heterotremus AB325521 IND 
rl- P. heterotremus AB325520 IND 
— P. heterotremus AB325518 IND 
P. heterotremus AB325517 IND 
P. heterotremus AB325519 IND 
PphC01ST5 KC859942 
PphC01ST7 KC859944 
PphC01ST8 KC859945 
PphC01ST6 KC859943 
PphC01ST3 KC8B9940 
PphC01ST1 KC859938 
,PphC01ST2 KC859939 
PphC01ST4KC859941 
PphCOiST9 KC859946 

P. pseudoheterotremus EF446313 THA 

- P. pseudoheterotremus EF44631 7 THA 
P. pseudoheterotremus EF446316 THA 
P. pseudoheterotremus EF014339 THA 
P. pseudoheterotremus EF446314 THA 

- P pseudoheterotremus EF44631 5 THA 
P. pseudoheterotremus JQ79681 5 THA 

87 1 1 PphCOl MW5 KC859936 
76lPphC01MW6 KC859937 
— P westermem ; AB354225 THA 
P wes(ermani AF333281 KOR 



93 



41 



66 



98 



71 



99 



I P. westermani U97207 JPN 
r P wesfemrara AY1 40686 PHL 
~96l P westermani U97213 PHL 



Fig. 3. Maximum likelihood tree based on partial cytochrome c oxi- 
dase subunit I (cox7) gene sequences. Sequences of P. heterotre- 
mus, P. pseudoheterotremus, and P. westermani (outgroup) ob- 
tained from GenBank are indicated with accession number and 
country code (ISO 3166-1 alpha-3 codes). Paragonimus sequenc- 
es of this study are presented in bold (KC859926, KC859936, 
KC859937, KC859931 , KC859927, KC859933, KC859935, 
KC859938-KC859946). The sequences were deposited in Gen- 
Bank numbers as shown in Table 1 . 



ent localities in Thailand and Myanmar were almost complete- 
ly (98-99%) identical. Phytogenies inferred from all available 
partial coxl sequences contained several well-supported clus- 
ters (Fig. 3). Sequences from Indian P. heterotremus formed a 
sister group to sequences from P. pseudoheterotremus-tfpe meta- 
cercariae. Sequences of P. heterotremus from Thailand, Vietnam, 
and China formed a separate distinct clade. A metacercaria 
from Phitsanulok Province was distinct from all others. Table 
2 shows the ranges and means of percentage differences be- 
tween these clusters, and that P. heterotremus from northeastern 
India are closest to sequences from metacercariae of the P. 
pseudohetervtremus-tfpe from Thailand. 

DISCUSSION 

In this study, we explored the molecular relationships of 
metacercariae conesponding morphologically with those de- 
scribed for P. heterotremus and P. pseudoheterotremus. Samples 
used were collected in different parts of Thailand including the 
Thai-Myanmar border area. Metacercariae corresponding to 
the P. heterotremus type were found in northern and central parts 
of Thailand and Myawaddy Province, Myanmar. Those of the 
P. pseudoheterotremus type were found in Surat Thani Province, 
Thailand and in Myawaddy Province, Myanmar. Members of 
this complex are clearly widespread in Thailand and the Thai- 
Myanmar border region, and this must alert us to possible hu- 
man infections in those areas. The ITS2 sequences were near- 
identical among all metacercariae sampled, and with sequenc- 
es available from GenBank. This reflects the situation seen in 
the P. skrjabini complex, which exhibits far less variation in ITS2 
sequences than in mitochondrial genes [16]. 

Although P. pseudoheterotremus was proposed as a species 
distinct from P. heterotremus, the phylogeny in Fig. 3, based on 
mitochondrial coxl sequences, shows that the true situation is 
complicated. P. heterotremus from northeastern India is sister to 
Thai metacercariae of the P. pseudoheterotremus type, and 1 meta- 
cercaria from Phitsanulok Province in Thailand is distinct from 
all others. The considerable variation in coxl sequences found 
in Thailand indicates that further work is required to demon- 
strate the full diversity of P. heterotremus and related forms in 
Thailand. 
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